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Abstract. We use the technique of electro-magnetically induced transparency (BIT) to probe 
the effect of attractive dipole-dipole interactions in a highly excited Rydberg gas. The transient 
character of the BIT response is investigated by rapidly scanning the probe laser through 
resonance. We characterize the resulting cooperative optical non-linearity as a function of 
probe strength, density and scan direction. For the 58D5/2 Rydberg state, an atom density 
of 1.6 X 10^° cm^^ and a positive frequency scan we measure a third-order non-linearity of 
10 ^ m^V ^. For the reverse scan we observe a second order non-linearity of 
X*^) = 5 X 10 ^ mV ^. The contrasting behaviour can be explained in terms of motional 
effects and resonant excitation of Rydberg pairs. 
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1. Introduction 

The excitation of laser cooled atoms to highly excited Rydberg states has opened many 
new research directions in quantum information [Q], many-body quantum physics |i2| and 
plasma physics (31. The use of electromagnetically induced transparency (EIT) to probe 
cold Rydberg gases flU provides a complementary technique to ionization for studying the 
properties of strongly interacting Rydberg gases [j5] and opens new perspectives in non-linear 
optics jU. The effect of the strong dipole-dipole interactions between Rydberg atoms on EIT 
is to produce a non-local optical non-linearity El. As such non-linearities can be large they 
are of considerable interest in the context of non-linear optics at the single photon level [Ulllll. 

Electromagnetically induced transparency arises from interference between an atom and 
two coherent optical fields. A three-level atom with ground state l^f), intermediate excited state 
|e) and Rydberg state |r), shown schematically in figured] (a), is driven by a strong coupling 
laser and a weak probe laser with Rabi frequencies Hp, He respectively. On the two-photon 
resonance, the atom evolves into a dark state, 

which no longer couples to the probe field, leading to transparency. For a pair of atoms 
separated by a distance R, dipole-dipole interactions between the Rydberg states cause a 
shift of the Rydberg pair-state energy level. For V{R) > 7eit> where 7eit is the bandwidth 
of the EIT resonance, this leads to a blockade of excitation of the Rydberg pair state [[TOl . 
The blockade mechanism modifies the non-interacting two-atom dark state |0) as f2p is 
increased, allowing population of |e) which resonant couples to the probe laser [|6] HI]. The 
result is a suppression of transmission on the two-photon resonance, and hence a cooperative 
optical non-linearity mediated by interactions with surrounding atoms. This non-linearity is 
enhanced as more atoms are included in the blockade region which leads to a quadratic density 
dependence 0. 

The optical response of Rydberg ensembles is sensitive to the nature of the dipole-dipole 
interactions between Rydberg atoms, which can be either repulsive or attractive and depend on 
the alignment angle between the dipoles. Rydberg interactions in cold gases has been the topic 
of in depth reviews, both experimental [[T2l and theoretical dH. In previous work [|6l the focus 
was on riSi/2 Rydberg states with zero orbital angular momentum, in the long range limit 
where the dipole dipole interaction is described by a repulsive, spherically symmetric van der 
Waals type interaction V{R) = —Cq/R^. While this case allows a direct interpretation of the 
underlying physics (ruling out alternative mechanisms such as ionization), larger interactions 
and consequently larger optical non-linearities may be accessible by employing Rydberg 
states with non-zero angular momentum. In this paper we focus on the 58D5/2 Rydberg 
state, coupling to the M = mji + mj2 = 5 pair state which experiences a strong attractive 
interaction with a van der Waals interaction strength of Cq ^ 150 GHz /im^ llT3l compared 
to Cq ^ —140 GHz /im^ [[T4| for the nearby 60S'i/2 state. A consequence of the attractive 
interactions is to introduce a temporal dependence into the EIT due to motional effects in the 
atomic cloud, which is explored through studying the spectra as a function of both frequency 
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Figure 1. (a) Three-level atom system. The probe laser detuning Ap = Wp — LOeg is scanned 
across the \g) to |e) transition whilst keeping the coupling laser locked on resonance with |e) to 
|r). (b) Dipole-dipole interaction potential for the 58£'5/2 pair state. The dashed line indicates 
the blockade radius _Rb below which only a single atom can be excited to the Rydberg state on 
resonance. States with an energy shift of —3 to —5 MHz are highlighted, corresponding to an 
interatomic separation ofi? = 5.6 — 6.1 jiva. 



and time. 



2. Experiment 

To determine the optical non-linearity for the case of attractive interactions, we measure 
transmission through a laser cooled ensemble of *^Rb atoms using the experimental setup 
described in ref. [4]. Atoms are loaded into a magneto-optical trap (MOT), reaching a peak 
density of p = 1.6 x 10^° cm"'^. A 10 ms optical molasses reduces the temperature to 
20 fiK, after which the atoms are prepared in the 5s {F = 2,mp = 2) state (\g)) 

by a 3 ms optical pumping pulse on the 5s F = 2 — )■ 5p ^P3/2 F' = 2 transition with 
polarisation. The density of the atomic sample is controlled by turning off the repump light 
resonant with the 5s F = 1 — )• 5p ^P3/2 F' = 2 transition during the optical pumping 
pulse, causing atoms to collect in the F = 1 hyperfine ground state. This allows the density of 
atoms in \g) to be reduced without changing the cloud size. Following a 1 ms free expansion, 
EIT spectroscopy is performed using counter-propagating probe and coupling lasers, both a^- 
polarized. The coupling laser is a frequency doubled diode laser system at 480 nm, which is 
stabilized to the 5p ^P3/2 -F' = 3 — )■ 58d ^05/2 transition using EIT in a room temperature 
Rb cell [[Bl [I6l. This is focused to a l/e^ waist of 215 ± 10 /im with a power of 80 mW, 
giving a peak Rabi frequency of nc/2n = 4.6 MHz. The probe laser is scanned across the 5s 
^Si/2 F = 2 (\g)) — 7- 5p ^P3/2 F' = 3 (|e)) transition over a period of 960 /is. The detuning is 
linearly ramped either from Ap/2n = —20 — )■ +20 — )■ —20 MHz corresponding to negative- 
positive-negative (NPN) or +20 —7- —20 — +20 MHz corresponding to positive-negative- 
positive (PNP) frequency order This double scan technique provides useful information on 
both atom loss H and hysteresis in the experiment. The probe beam is focused to a waist of 
160+10 yum, with data recorded for powers varying from 5-40 nW corresponding to peak Rabi 
frequencies in the range ^lp/2ir = 0.3-1 MHz. Transmission through the cloud is recorded 
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using a fast photodiode with a 20 MHz bandwidth. The optical path length £ is determined 
from time of flight imaging of the cloud ifFTl . giving i = 0.9 ±0.1 mm. 

3. Temporal dependence of the EIT lineshape 

Theoretical treatments of the optical non-linearity arising from dipole-dipole interactions 
in EIT lITl [m have so far assumed the sample to be frozen, with no modification of the 
interatomic separations in the atom cloud due to interactions. However, even at 20 ixK the 
motional effects from the attractive interactions can be significant on timescales of tens of 
microseconds [[T9l l20l . Figure [11(b) shows the dipole-dipole interaction potential V{R) for 
the 58L'5/2 pair state as a function of interatomic separation. For the experiments presented 
below, 7EiT/2vr ~ 1 MHz corresponding to a blockade radius of i?b = ^/Cq/jevt — 7 fim 
which is indicated by a dashed line. Comparing this to the average interatomic separation 
given by (R) = {5/9)p^^^^ EB, for a density of 10^° cm'^ (R) ~ 2.5 /im so the blockade 
effect is expected to be clearly observable for the experiment parameters. For atoms initially 
separated by 7 /im, the collision time is approximately 20 jj,s [l20|. resulting in ionisation 
of at least one of the two atoms. These two effects of motion and ionisation introduce a 
temporal dependence into the EIT spectra, making the optical non-linearity a function of not 
only optical fields, but also the timescale over which they are probed. 

In addition to the resonant blockade behaviour, the EIT spectra are also strongly 
influenced by the off-resonant excitation of the anti-blockaded pair states at Ap = V{R)/2 
that have a separation R < i?b, and therefore ionise more rapidly. As an example, states with 
V{R)/27r = -3 to -5 MHz (corresponding to a detuning of Ap/27r = -1.5 to -2.5 MHz) 
are highlighted in figure [11(b), with interatomic separations of i? ~ 5.8 /im, with a collision 
time of around 8 /is. This has two effects, firstly creating ions off-resonance that can then 
affect the behaviour on the EIT resonance, and secondly it reduces the number of short-range 
pairs in the medium. Off-resonant excitation of these anti-blockade states has previously been 
used as a method of mapping out the distribution of nearest-neighbours within the atom cloud 

m. 

To demonstrate these temporal effects, EIT spectra are recorded at a density of p = 
0.9 X 10^° cm^'^ with f2p/27r = 0.9 MHz for a positive scan across the resonance for a range 
of scan speeds, shown in figure |2](a). For the slowest frequency ramp of 50 MHz/ms, the laser 
scans across the EIT resonance in 25 fis. Around Ap/27r = —2 MHz, corresponding to the 
excitation of the anti-blockaded states highlighted in figure [I](b), the transmission increases 
giving a broadened and very asymmetric EIT resonance which has a narrow absorption profile, 
with atoms appearing to be lost rapidly around Ap/2n ~ 1 MHz. This spectrum is consistent 
with significant ionisation of the off -resonant pair states, which appears to dominate over 
the EIT suppression expected due to blockade, seen from a broadening of the EIT profile. 
Figure [2] (b) shows calculations of the pair-state population using the pair model detailed in 
ref. as a function of interaction strength for r2c/27r = 2.4 MHz to match the experiment 
parameters. This demonstrates the suppression of the resonant pair excitation due to blockade, 
whilst highlighting importance of the off-resonant excitation of the anti-blockade states for 
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Figure 2. Temporal effects in EIT. (a) Measured EIT spectra recorded at 30 nW (rip/27r = 
0.9 MHz) for a range of scan speeds. For the slowest scan (50 MHz/ms) there is evidence 
of broadening and loss consistent with ionization or motional effects starting at the detuning 
corresponding to pair states highlighted in figure [T] (b), whilst at faster speeds the spectra 
appear consistent. Fitting the high-scan rate data gives VIc/^-k = 2.4 MHz and a relative 
two-photon linewidth of 200 kHz. (b) Rydberg pair population calculated using the pair- 
model described in ref. (|6l calculated for different interaction strengths. This shows resonant 
excitation of the pairs with R < before scanning across resonance, explaining the 
ionisation in (a). 



negative detunings. For large interactions, the off-resonant pair excitation is suppressed due 
to being far-detuned, which is why the effects of the pair states only become observable 
around Ap/27r ~ —2 MHz despite the average interatomic separation giving a shift of several 
hundred MHz. From the groundstate density, the probability of finding atoms separated by 
5.8 /im is around 0. 1 %, corresponding to a pair state density around 10^ cm~'^ at this detuning, 
which can then ionise and broaden the resonance. As the scan speed is increased, the spectra 
in figure |2] (a) become symmetric with approximately constant transmission on resonance, 
showing that loss due to ionization is suppressed as there is not enough time for a significant 
number of ions to be created during the scan across resonance. All data presented below is 
recorded for a scan speed of 80 MHz/ms, fast enough that there is no evidence of loss or 
asymmetry across the first EIT resonance but slow enough to enable these dynamic effects to 
be relevant. 

Despite there being no observable loss or asymmetry across the first resonance at this 
intermediate scan speed, there may be a small ion fraction in the cloud. For S-states as studied 
in [IHl, it is easy to assess the role of ionization in the suppression data as this leads to a red- 
shift of the EIT resonance. For the /^-states however, the presence of an ion in the sample 
breaks the quantisation along the z-axis and will project the Rydberg atom into one of the 
\mj\ = 1/2, 3/2 or 5/2 states which have scalar polarizabilities of = —137, 111, 607 
MHz/(V/cm)^ respectively, calculated using the Numerov method [[23]|. To estimate the effect 
this mixture of polarizabilities has on the lineshape, a model was developed to calculate the 
lineshape for a uniform density of atoms, from which a fixed fraction were chosen as ions. For 
each of the remaining atoms, the total field electric field due to the surrounding ions is used to 
find the Stark-shift assuming a random choice of \mj\, and the total lineshape calculated by 
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Figure 3. Role of ions in EIT suppression, (a) Simulated Hneshape as a function of ion 
fraction for a random atomic distribution. Due to the sign difference in ao of the different 
rrij states, the lineshape isn't simply red-shifted as expected for S1/2 states but instead the 
resonance disappears, becoming asymmetric in the wings, (b) Transmission data recorded 
at peak density showing strong suppression of the EIT for a strong probe as expected from 
interactions, however this shows good agreement with the ion model for a 3 % yield. 



summing over the susceptibility of each atom. 

The results of the ion model are plotted in figure |3](a) as a function of ion density within 
the cloud. This shows that even for a very weak ion yield, the Stark-shift of the two-photon 
resonance significantly modifies the resonant transmission. This doesn't simply correspond 
to a red-shift, instead the line is split due to the varying rrij components, reducing the 
transmission whilst becoming increasingly asymmetry either side of the original resonance. 
This is contrasted in figure |3] (b) by experiments performed for a positive frequency scan at 
the peak density of p = 1.6 x 10^° cm~^ where ionization effects should be maximized. 
For the weak probe data the EIT lineshape is unperturbed, whilst at np/2n = 0.7 MHz 
the transmission feature is strongly suppressed as expected due to interactions. However, it 
looks similar to the Stark detuned profile expected for a 3 % ion fraction. One caveat to 
this comparison is that the model assumes ions are present at all times in the scan, and not 
created dynamically. As seen from figure |2](b) the off-resonant pair-states, and hence ions, 
are most likely to be excited close to the EIT resonance. In addition, the blockade mechanism 
suppresses the excitation of close Rydberg pairs which should prevent ion generation during 
the scan. As there is no way of measuring the ion yield in the current experiment, it is 
difficult to assess the relative importance of ionization over blockade for the positive scan 
data, however it is still interesting to consider a system that gives a non-linear optical response 
due to other mechanisms such as ionization. This could be extended to the idea of performing 
EIT on highly excited Rydberg states and using the Stark-shift to create an ion-blockade 
[|24l which would switch the system from transmission to absorption, making a sensitive ion 
detector. 
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4. Frequency dependence of EIT suppression 

To analyse the dynamical response further we compare the behaviour of a positive frequency 
scan to the negative frequency scan. Figure |4] shows data recorded for both NPN (a) and 
(b) and PNP (c) and (d) scan directions for weak r2p/27r = 0.3 MHz (a) and (c) and strong 
0.9 MHz probe fields (b) and (d) at the peak density of p = 1.6 x 10^° cm^'^. In the weak- 
probe regime the dark- state |0) = \g) and the interaction effects play no role in the behaviour 
of the system. Comparing data in (a) and (c), the scans show symmetric EIT traces with 
no hysteresis between the first and second scan, as expected. This confirms the ground-state 
density is not varying over the scan duration. For the strong probe data for the NPN scan 
direction in (b), the first resonance is strongly suppressed as observed above, however the 
reverse scan in the negative direction shows almost complete recovery of the EIT comparable 
to the weak-probe trace in (a), with a slight reduction in the optical depth relative to the first 




Figure 4. Rydberg EIT spectra for the 55*1/2 ^P3/2 58Z55/2 system for NPN (a) and 
(b), and PNP (c) and (d) scan directions at probe fields of rip/27r = 0.3 MHz (a) and (c) 
and 0.9 MHz (b) and (d). For the weak probe data (a) and (c) symmetric EIT is observed on 
the two photon resonances. At higher probe Rabi frequency (b) and (d), there is evidence 
of suppression of the EIT peak and asymmetry dependent upon the direction of the scan 
through resonance. For a positive scan there is a strong suppression, which is revived in the 
reverse scan at later time (b). The negative scan in (d) reveals loss in the second scan. For all 
traces, the frequency corresponding to excitation of the anti -blockade states from figure [T](b) 
is highlighted. 
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scan. Contrasting this to the opposite PNP case (d), on the initial scan the transmission is 
suppressed but only by around half that observed in the NPN data. The second scan in the 
positive frequency direction however shows significantly greater suppression of the resonance, 
in addition to around a 10 % change in the background absorption from loss of atoms in \g). 
As well as the line- shape, the frequency of the two-photon resonance was compared for 20 
datasets for a variety of probe powers. This shows no systematic shift for the strong probe 
data over the ±300 kHz noise observed in the weak probe limit arising from fluctuations in 
the frequency of the coupling laser. 

This data reveals the anti-blockade states play an important role in determining the 
magnitude of the EIT suppression, with the short-range pairs identified in figure \T\ (b) 
highlighted here in grey for each scan. This clearly shows that the EIT suppression for 
the strong probe data of (b) and (d) is significantly larger when scanning from negative 
detuning across the short-range pair states compared to the suppression scanning from positive 
detuning. The enhancement is most likely due to ionisation of the short range pairs causing an 
irreversible ion-blockade to dominate over the coherent excitation blockade due to the longer 
range of the Coulomb potential. The pair-state ionisation also explains the hysteresis in the 
datasets, as for (b) the NPN scan first creates a strong suppression by excitation and ionisation 
of the short-range pair states. This reduces the number of atoms with separations of R < i?b, 
leading to the observed recovery in the reverse scan 0.5 ms later. This hysteresis in the 
strong probe data is interesting, as the optical response is now dependent upon the experiment 
sequence modifying the atomic distribution. This is similar to direct observations of ion-yield 
for attractive states excited with either a red or blue detuning relative to the Rydberg energy 

mm. 

5. Optical Non-linearity 

As seen from the data in figure |4]the suppression of the first scan for the PNP data is due 
to the dipole-blockade mechanism, whilst for the NPN scan it appears to be enhanced by 
ionisation of the short-range pair states. To characterise the signatures associated with these 
two regimes, the magnitude and density dependence of the resonant optical non-linearity is 
measured. Data is taken for a range of probe powers and densities for both scan directions. 
For each set of parameters, a total of 20 traces are used to determine the transmission T on 
the first EIT resonance. This is converted to the complex part of the susceptibility at the probe 
frequency using the relation Im{x} = —\og^{T)/k£, where k = 27r/Ap is the wavevectorof 
the probe laser. 

The measured susceptibility is plotted as a function of the probe electric field £ in 
figure |5l for both the NPN (a) and PNP (b) scan directions. In the high-density data in 
(a) the non-linearity is seen to saturate around 20 V/m (ilp/27r = 0.8 MHz), after which 
it reverses direction. This turning point appears due to loss of atoms as the laser scans across 
the resonance, and places a limit on the largest probe power that can be used. All of the curves 
in (a) scale consistent with a third-order non-linearity. The function x = X^^^ + x^'^^£'^ is fit 
to all the datasets, giving a peak value of ~ 5 x 10"'' m^V~^ for the highest density data. 
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This represents a very large non-linearity for an atomic sample, comparable in magnitude 
to that observed in slow light experiments at BEC densities [|26ll . For the PNP data in (b), 
the maximum suppression of the resonant transmission is lower, and there is no evidence 
of saturation even at strong fields and high density. The susceptibility reveals a linear scaling 
with E consistent with a second order response and is fit to x = x^^^+x^^^f^, giving a peak value 
of x*^^-* = 5 X 10~^ mV^^. This difference in the non-linear scaling for the two scan directions 
is consistent with there being different mechanisms (blockade or ionisation) responsible for 
the suppression for each scan direction, as discussed above. 

Comparing this data to the theoretical results obtained from the many-body Monte-Carlo 
model of Ates et al. the model predicts a third-order non-linearity which saturates as f2p 
approaches fic for the suppression due to blockade from repulsive interactions. However, 
from the data a third order non-linearity is observed for the case of enhancement due to 
ionisation, whilst the suppression due to dipole blockade gives a second order scaling. This 
discrepancy requires a more detailed study involving direct measurement of the ion fraction 
for each scan direction to explore the importance of these temporal effects compared to the 
theory that considers a static atomic cloud to try and elucidate the reason for this difference 
in the order of the optical non-linearity. Despite the data in (a) and (b) showing different 
non-linear scalings, the maximum values of the complex part of the susceptibility for a given 
density are comparable. This is limited by two factors - the first is the absolute optical depth 
of the sample in the absence of the coupling laser, as this sets the maximum value of x in 
the experiment, and the second is the number of atoms per blockade sphere, J\f. As more 
atoms are added to the blockade, the dark- state fraction reduces proportional to l/J\f, leading 
to a non-linear density dependence. The density dependence of the non-linear susceptiblity is 
shown in (c) and (d) for the NPN and PNP data respectively. This reveals a clear quadratic 
scaling of x^^^ with density and an initial quadratic dependence for x^^-* that saturates at higher 
densities. These observations are consistent with the expected density scaling found in ref. 0. 

6. Conclusions and Outlook 

In summary, we have shown that the optical response observed for EIT in an attractive gas is 
strongly dependent upon the hysteresis in the system due to motional effects and ionisation 
of the red-detuned anti-blockaded states. For a positive frequency sweep, the resulting third- 
order non-linearity is very large due to the enhancement from the ionisation of the short-range 
pair states. A negative frequency sweep however reveals a second-order non-linearity due to 
the blockade effect. 

In future work, the attractive potential will be characterized in the dispersive regime 
for a single blockade volume to explore its application to single-photon quantum gates. An 
additional parameter accessible for non-zero angular momentum states is the dependence on 
the alignment of the dipoles [|27l . This could be explored with the application of a weak 
electric field, potentially enabling an external switch to turn the blockade-mechanism on or 
off. 
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Figure 5. Non-linear susceptibility. Imaginary part of the susceptibility as a function of the 
probe electric field 8 for the positive (a) and negative (b) scan directions for p = 0.4 (O) 0.7 
(□) 1.0 (•) 1.6 (■) xlO^" cm^^. Data in (a) and (b) is fitted to third- and second-order non- 
linearities respectively, and the resulting density dependence plotted in (c) and (d). Both x^^^ 
and x^'^'' display a quadratic density scaling, consistent with pair-wise interactions. 
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